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CHAPTER 1 
INTRODUCTION 


In the past two decades, the measurement of S (Scattering) parameters has 
become progressively more important. This importance stems both from the ease of 
the basic measurement, particularly over a broad frequency range, and the usefulness 
of this parameter set in the design process. Measurements have been made to below 
1 MHz and are slowly being extended into the millimeter wave spectrum (though 
measurements in this range are almost by necessity done in an S— parameter format). 
We, not unsurprisingly, find the calibration process and standards at the heart of 
such a measurement system. 

To calibrate such a system for reflection and transmission measurements, 
Hewlett-Packard Company (HP) provides a set of standards (originally from Maury 
Microwave) and related data for calibration. Using the HP procedures, several 
researchers have found anomalies in the resultant data. A not uncommon phenomena 
is an air— line with a reflection that shows gain. v Admittedly, this gain exceeds unity 
by less than one percent. However, somewhat by accident, researchers have used 
slightly different calibration data to represent the standards and found substantially 
improved results for the airlines. 

The purpose of this report is to investigate the full calibration process for the 
existing HP measurement system. We will begin by reviewing the basic calibration 

process and measurement system. The approach taken by Hewlett-Packard will be 
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discussed as the methods are developed. A full model for the system will also be 
developed. This model may be used to develop full system calibration, including the 
calibration of the HP "open" circuit reference. This approach has the additional 
advantage of not requiring the "open" for the calibration, a fact critical to certain 
measurement environments such as the measurement of dielectric and conductive 
material samples. In such instances, the calibration must often be done using some 
alternate procedure such as that presented in this report. More importantly, the 
basic formalism of the modeling enables a user to develop a calibration procedure 
tailored to a specific need where the HP procedure is inadequate. Our premise is to 
improve the HP calibration data and then use the fundamental HP procedure when 
possible. This will provide a high degree of accuracy if the correct calibration data 
is available, usually in excess of the accuracy specifications of the HP 8510 network 
analyzer system. 

Sample results will be presented which validate the method and demonstrate 
the potential errors using HP's process and data. This will include the effects of 
both fixed terminations and sliding loads on the choice of calibration data used to 
represent the HP "open" standard. Limitations of the approach will be discussed in 
addition to assumptions built into the calibration and accuracy process. Associated 
effects of the assumptions about switching networks have led to alternate procedures 
requiring a system modification (Speciale, 1977b). Since this latter work is not 
directly related to the report objective, no consideration will be given to this 
alternate procedure at this time. 
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CHAPTER 2 

DEVELOPMENT OF 8-PORT MODEL 
AND DESCRIBING EQUATIONS 


The measurement of the S parameters of a device under test (DUT) is to be 
performed using a HP network analyzer system. The development to be presented is 
applicable to both the HP8410 and HP8510 series systems, excluding the new TRL 
calibration technique. The basic variations in the measurements are due in part to 
each of the following in addition to some other minor factors: 

1. Imperfections in coaxial switches used in the system, particularly 

causing variations between switching periods; 

2. System noise; 

3. Generator noise; 

4. Frequency stability (improved with phase lock); 

5. Nonlinearity in detectors (particularly phase detection near the 

limits); 

6. Signal leveling; 

7. Quantization errors (in A-to— D conversion). 

The basic system is shown in Fig. 1 as used for measurement of s and s^. 
To measure s^a and s^, either the DUT may be reversed or the source may be 
switched to the right end of the system. The latter is used for the HP systems and 

is a random source of error, which may be avoided only with a system modification. 
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The modification (Speciale, 1977a) will not be considered in this report. We will 
find later that the calibration may be done for the system as shown and repeated 
independently for the system with the source switched to the output. The basic 
procedures for the two sets of measurements do not interact. 

Though Fig. 1 represents the fundamental layout, the actual measurement 
process is complicated somewhat by possible crosstalk between channels of the 
network analyzer. A complete 8— port model of the 2— port measurement system with 
the DUT connected is shown in flow-graph form in Fig. 2. If we observe the details 
of the measurement system in Fig. 1, we find 8 ports to the outside world: the two 
source connections, the DUT connections, and the 4 network analyzer terms. The 
latter have been referred to as terms since they may actually represent voltages 
internal to the network analyzer which are processed by the network analyzer for the 
S— parameter estimation. 

The network of Fig. 2 may be substantially simplified using either of two 

basic observations. Since the terminations on the network analyzer ports and source 

ports (neglecting switching errors) are to be constant at a given frequency (not 

necessarily zero), they may be imbedded in the S— parameters of a simpler system 

description. This is a common practice when an n— port is to be used as an 

(n— m)— port with the extra m ports terminated in fixed loads. The computation 

requires a bit of matrix algebra. To demonstrate this effect, let us assume the 

parameters of the 8— port are described by s'- ■ and port k is terminated in I\. Then 

ij a ; 

the modified ■ are given by 




s ik s 'kj^k 
1 ‘ s kk T k 


(la) 
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and 


(lb) 

(lc) 

These formulas may be used for each port k which has a constant termination until 
the full system has been simplified. 

An alternative is to note that a similar result may be obtained by using the 
general form of S— parameters and defining the normalization impedance at each of 
the constant termination ports to be the constant termination. This approach results 
in the same form. The actual numbers obtained from the original are unimportant 
since the true numbers are not known. In fact, the desired result is the form 
obtained and not the actual numerical relationship of the final model numbers to the 
original. This form for the 8-port system which describes the HP network analyzer 
systems is given by the simplified flow graph of Fig. 3. For this form, the final 
system S-parameters will be denoted by sj ., and the relationship to Fig. 2 will be 
ignored. In this simplified system, only 24 parameters are nonzero. Indeed, if all of 
the parameters existed, they would have been found to combine in the calibration 
process to result in another set of 24 terms. 

Before we consider the full calibration problem, we will consider the 
fundamental problem of reflection coefficient measurement of a one— port. In this 

case only s ^ of the DUT is nonzero. We will denote this value simply by s. For 
this measurement, we may set (the second source) to zero and consider only the 
ratio r = (t j / i j). In an ideal system, T would be equal to s. Our job is to 


S kj 


Ik. 


s k T k 


or 


sS ik 


1 - 


s kk T k 
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determine the T and s relationship so that 5 may be computed from a measured 
value of T. 


In terms of the system S— parameters and s, the measured T displayed by the 
network analyzer will have the value 


or equivalently 


r = 


5 s 4- fs s s s 

*31 _ J 1 


’n 31 


s s 01 ) s 


$ 6) 1 “f” ( $ ! & 1 1 S 


5 + (s s t - S 


’ 21 


r = 


*24 41 


Oq + dj s 
1 + s 


s 

44 


S 21 ) 5 


( 2 ) 


where the • are constants and the numerator and denominator have been divided 

ij 

Q 

by s 2J which represents the incident coupling to the incident sample port. Since 
only the ratio T is desired, we may simply normalize all signal paths departing the 
source by the level of coupling to the incident sample port if desired and reduce the 
number of unknown coefficients by one. 

As already suggested, we need to determine the coefficients of (2) in order to 
compute s from the measured T. If we assume that ideal terminations are available, 
then clq is the T obtained with a perfect terminating load attached. For the fs of a 
short and open given by T and T respectively, then 

o U 


Q 1 = 


^0 + V “0- 2T c T s 


and 


2 a n - r - r 

1 0 0 s 

1 ~ r - r • 

0 s 

With these terms, then s is simply computed as 

a 0 ~ r 

5 ~ -a j + T [ T • 


(3) 
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Thus the basic process for a reflection coefficient s is rather straightforward. 
The problem arises when we realize that the only near-perfect termination is the 
short circuit. The load and open are essentially impossible to construct as perfect 
terminations. This means that a procedure must either be developed which pursues 
another process of calibration, or the load and open terminations must be made 
known in some manner. These processes will be the discussion of latter chapters of 
this report. We stop at this point of the simple reflection problem and begin the 
more general problem and the necessary calibration procedures in the next chapter. 

The full development- of the next chapter requires to use of the basic equations 
for the model of Figure 2. These equations are given by 


l l - 

S 21 

b sl 

+ 

c 5 

s 24 

b i 

+ 

s 25 

b s 

+ 

s 28 

b s2 

(4a) 

h = 

s 3l 

41 

+ 

s s 

s 34 

b , 

+ 

s 35 

»* 

+ 

c S 

8 38 

b s2 

(4b) 

l 2 = 

s s 

8 71 


+ 

s s 

s 74 

b , 

+ 

s 75 

b 2 

+ 

c 5 

s 78 

b sS 

(4c) 

*£ = 

S 61 

b s, 

+ 

s s 

s 64 

b i 

+ 

sS 65 

»* 

+ 

s s 

s 68 

b s2 

(W 

a l = 

S 41 

b sl 

+ 

S tt 

b , 

+ 

4s 

b l 

+ 

s 48 

b s2 

(4c) 

a 2 = 

4, 

b sl 

+ 


b , 

+ 

4s 

4 

+ 

s 58 

b s2 

(4f) 



b l 

— 

s n 

a i 

+ 

s 12 





(49) 



b 2 

= 

S 21 ' 

•t 

+ 

s 22 1 





(4b) 


These equations are combined in the next chapter to obtain a complete description of 
the measured S-parameters. 
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CHAPTER 3 

CALIBRATION USING IDEAL STANDARDS 


As described in the previous chapter, a rigorous development of the system 
model and resultant equations is easily done. We shall complete the equation 
development in this chapter and begin the consideration of full system calibration 
with the use of perfect calibration standards. This development will be extended in 
later chapters. 


3.1 Describing Equations 

As shown in the previous chapter, the simplified eight-port model of the 
network analyzer has eight describing equations. These equations may be combined 
to obtain equations for only the four measured terms i j, t ^ and It is 


t 


convenient to define (— 


'1 


b s2=° 


), h 


b s2 =0 


), (i 


‘2 


b sf° 


), and (- 


t, 


l 2 


b sl=° 


) , denoted 


s® i , s'21 , s"^) and s' 22 respectively (the m denoting measured) to represent the four 
ratios are measured directly by the network analyzer, to obtain 


m 


m 


m 


9 11 


- { 


S s 0i + (s s 


31 

s s 


s s a s 
34 s 41 s 31 s 4i 


) 8 11 + (s 


11 

s s 


34 s 51 s 31°54 
s s 


) « 


12 


+ ( s 35 s 41 s 31 s 45 S 21 + ( s 35 s 51 s 31 s 55> 5 


22 


, r S / S 3 S S \ , S / S S S S \ 

+ 1 s 34 51 s 45~ s 41 S 55^ + s 35 ( s 4t s 5f s 51 s 4$ 


, "S / s s s s 
+ S 31 ( s 44 s 55 s 45 s 54 


)] 4 S ) / {< 


.6 , / S 0 S S \ 

21 + ( s 24 s 4i~ s 2l s 44 ) 5 


11 
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+ ( s 24 s 51 s 2 I s 54 * s 12 + (' s i'r' s // _5 ' S “ 5 


r) s ( 


25 41 21*45' 9 21 

+ ( sS 25 s 5r sS 21 s 55) s 22 + K 24 ( s 51 s 45~ s 41 s 55) 

+ s 25 ( 3 41 8 54- 3 51 s 44) + 4l (5 ^ s 55“^5 s 5^^ A s } 


(5a) 


s m 

S 21 


* 61 


+ {si 


64 s 41 S 6l s 44> s ll 


) s ]{ + 5 


and 


64*51 s 61 s 54^ 3 i 
) s t 


+ ( s te s 4r4l s 4j) S 21 + ( s 65 s 5r s 61 s 55> S 22 
+ [ s 64 ^ S 51 sS 45~ S 41 S 55^ + 5 05 ( 3 41 3 54~ 3 51 3 44* 

+ S 61 (s ^ s 55 -5 ^5 5 J^J A s ) / {4/ + ( 3 24 sS 4r 3 21 sS 44) S 11 
+ ( s 24 s 5 l~ s 2 I s 54} s 12 + . ( s 254 f r s 21 sS 4 ^ S 21 


+ (-sirs'; 


~ 3 21 s 55^ 3 22 + I s 24 ( S 5 i S 45 s 4 I s 55) 


25 51*21* 55' 9 22 

+ S'tt r (s^ f Sr/ _ ir.S 


25 ^ 41*54 *51*44* + S 21 ^^ S 55 _S ^5 S 5^J A s } 


s 12 ~ { 5 5£ + 34 s 48~ s 38 s 44* S U + ^ S 34 s 58~ sS 38 sS 54* s 12 

+ ^ s 35 s 48~ s 38 s 45* S 21 + ( s 35 s 58~ s 38 s 55* s 22 
+ ^ ^5^5 -s ^55^ + s 35 

+ s 38 (5 ^ 5 lr^5 5 5^ A s } / { s 78 + ( 3 74 sS 4F 3 7& S 44* S 11 
+ ( s 74 s 5f s 78 3 54* s 12 + S 21 

+ ( s 7, 5? 8 58~ 3 78? a 55* 3 22 + l a 74 
+ s 75 ( s 48 s 54~ s 58 a 44* + s 78 A s } 


s m 

s 22 


{ 5 &? + ( s 64 s 48 3 68 s 44* S 11 + ( 3 64 3 58 s tfdT s 5^ s 12 
+ ( s 65 s 48~ s 68 s 45* S 21 + ^ s tf5 s 5^ _5 ^55^ 3 22 


(5b) 


(5c) 


12 


+ 

+ 

+ 

+ 

+ 


3 48 3 55^ + s 65 s 58 s 44^ 

s 68 ( s 44 sS 55~ s 45 S h^ A s } / { 

( 74 s 58 s 7 8 s 54’ s 12 + 75 s 48 3 78 s 4$ S 21 

( s s s s \ „ , r 5 t s s s s ^ 

( s 75 s 58~ s 78 s 55 ' s 22 + l 5 7^ ( 5 5<? 5 45 5 ^55^ 


5 v 

5 y , ) 5 


75 ^ ^°74°48 °78°44 } *11 
s s 


S / s s s 

3 75 ( s 48 s 54 s 58 s 44 


\ 1 5/55 

) + 5 7 R i S ll S K<~ S 


78 K 44 55 °45°54>' s 


)] 


(5d) 


The denominators and and numerators of these equations all have the form (a, + 
a 2 Sjj + a 3 Sj 2 + a A s 3I + a b s^g + a 6 A^). In addition, Eqs. (5a and b) have 
the same denominators as do Eqs. (5c and d). This last fact is useful in that the 
forward and reverse parameter sets may be calibrated and measured independently. 
In fact if the system user desires to perform the calibration only once, the first two 
equations may be used for the entire measurement process of a two-port by simply 
reversing the DUT in the measurement system. To summarize the equations for 
which the coefficients must be determined, we simply normalize the numerators and 
denominators of (5) by and as appropriate to obtain 


s m = 

a 0 

+ 

a i 

s n 

+ 

a 2 

s 12 

+ a 3 

S 21 

+ 


s 22 

+ 

a 5 


(7a) 

S 11 

1 

T 

A 

s n 

+ 

A 

3 12 

+ b 3 

S 21 

+ 

u 

s 22 

+ 

A 


si = 

c 0 

+ 

c i 

s n 

+ 

c 2 

s 12 

+ c s 

S 21 

+ 


s 22 

+ 

c 5 


(7b) 

S 21 

1 

T 

A 

s n 

+ 

A 

3 12 

+ b 3 

S 21 

+ 


s 22 

+ 

A 


.» = 

d 0 

+ 


s n 

+ 

d 2 

3 12 

+ d 3 

S 21 

+ 

L 

s 22 

+ 

d 5 


(7c) 

s 12 

i 

+ 

e i 

s n 

+ 

e 2 

s 12 

+ e 3 

S 21 
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e 4 

s 22 
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e 5 


s m - 

fo 

+ h 

s n 

+ f g 

3 12 

+ h 

S 21 

1A 

s 22 

+ f 5 


(7d) 

s 22 

1 

+ 

e i 

s n 
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e 2 

3 12 

+ e 3 

S 21 

+ 

e 4 

s 22 

+ 

e 5 
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In a perfect system, these four ratios are precisely equal to s 1} , S2j, s^, and ,s ^ 
respectively. Thus, in a perfect system coefficients a p c 3 , </^, and f ^ would be 
unity and all other coefficients would be zero. 

In a practical system, three nonideal effects appear. First, there may be both 
physical and electrical differences between the reference and test channels, resulting in 
values of a p c^, 62, and f ^ unequal to unity. Various user controls (e.g. amplitude 
and phase controls and line stretchers) are provided in a typical system to 
compensate for these differences, but the compensation can be made exact for only 
one frequency and one measured S parameter at a time. Thus, the use of 
swept-frequency measurements and S parameter test sets sacrifices accuracy for 
convenience. It is important to note, however, that accurate measurements require 

only that dp c^, d ^ , and be know as functions of frequency, not that they be 
set to unity a priori. 

The second non-ideal effect that may appear is the presence of constant offsets 
in the measurements (non-zero Oq, Cq, < 1 q, fg). These offsets may also be 
characterized during system calibration and their effects removed to produce accurate 
measurements. There are usually no user-accessible controls on the equipment to 
zero them out, although on a polar phase-magnitude (Smith Chart) display the 
centering adjustments effectively perform this function. 

The third non-ideal effect that may appear is crosstalk between parts of the 
network analyzer due to the finite directivity of the directional couplers and direct 
radiation (leakage) within the instrument. The mathematical result of this crosstalk 
is that any or all of the coefficients in the four describing equations may be 
non-zero. Once again, the values of these coefficients can be found during 
calibration; now, however, the four equations become coupled with respect to the 
unknown S parameters of the device, and numerical computation is required to find 
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the S parameters. (Further, the equations are nonlinear in the S parameters.) It 
must be borne in mind, however, that the effects of the crosstalk are small in a 
high-quality system, so that the a p c^, </g, an ^ fi terms are still the dominant 
terms in their respective equations. 

Calibrating the network analyzer system consists of computing the values of the 
coefficients in the describing equations, based on measurements of known devices. 
Possible approaches to doing this calibration are discussed in this chapter and the 
next, starting here with the simplest case that of using "perfect" standards. 

3.2 Sample Calibration Using Ideal Standards 

Assume that the following standard components are available to be measured by 
the network analyzer: 

1) Perfect shorts: two short-circuited transmission lines having reflection 

coefficients p = -1 + jO. 

2) Perfect opens: two open-circuited transmission lines having p = 1 + jO. 

3) Perfect loads: two terminated transmission lines having p = 0. 

4) Perfect through line: a transmission line having no loss, no discontinuities, and 
zero electrical length. 

5) Perfect half-wave through line: a transmission line having no loss, no 

discontinuities, and an electrical length of a half wavelength. 

The calibration procedure will be illustrated here by find the a-, b-, and c , 
i = [0,.. .,5], in the describing equations for and These two equations are 

decoupled from the remaining two describing equations with respect to the 
coefficients; hence, the two sets can be solved independently for the coefficients. 
(This statement is not immediately apparent, since the coefficients all depend on the 
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same set of S parameters {s^; 1 < ( i , j) < 8}. The proof will be by 
construction; i.e., by developing a procedure.) 

For reference, the two describing equations whose coefficients are to be found 
are repeated here: 


a* = 

a 0 

+ 

a l 

s n 

+ 

a 2 

s 12 

+ 

a 3 S 21 

+ 

“l 

s 22 

+ 

°5 4 S 

b ll 

T 

+ 

"*7 

s n 

+ 

~^2 

s 12 

+ 

b 3 S 21 

+ 


s 22 

+ 


3 " = 

c 0 

+ 

c i 

s u 

+ 

c 2 

s 12 

+ 

c 3 S 21 

+ 

u 

s 22 

+ 

c 5 A s 

S 21 

T~ 

+ 

Ti 

S 11 

+ 

~^2 

s 12 

T~ 

b 3 S 21 

+ 

U 

s 22 

+ 

; • 


(7a) 


(7b) 


3.2.1 Step 1 

Connect perfect loads to ports 1 and 2 so that 


Then 

and 



( 8 ) 

(9a) 

(9b) 


3.2.2 Step 2 

a) Connect a perfect open to port 1 and a perfect load to port 2. Then 


so that 


and 


S = 

*- 1 

(10) 

II 

<3 

“o + a l 

i + b j 

(11a) 

A 

ll 

c 0+ C 1 
1 + b j 

(lib) 
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so that 


and 


a perfect short to port 1 and 

a perfect load to port 2. Then 


S 

II 

1 

(12) 

s llb 

II 

a 

l 1 

(13a) 

s 21b 

c 0 C 1 
- r-^f 7 - 

(13b) 


After some algebra, the simultaneous solution of Eqs. (11a) and (13a) is found to be 
(recall that Qq is known) 

•t (An + s ub> - 1 An Ait 


•i = 


' 11a 


m 


(14a) 


'lib 


and 


A = 


2 a n - (s 


m 


m 


lla + s llb ' 


m m 

s m s l ib 


(14b) 


Similar expressions result for Cj and bj after solving (lib) and (13b) simultaneously. 
The two equations for calculating bj should give identical results, thus providing a 
check on measurement and computational techniques. Bear in mind, however, that 
an( ^ ( S 2 fb ) will b e ver y small compared with i s n a ) and (s^); 

consequently, values of b^ computed using (s^ fl ) and will be less corrupted 

by measurement noise. 


3.2.3 Step 3 

a) Connect a perfect load to port 1 and a perfect open to port 2. Then 

s - f 0 0] 

b - [0 1 J 


(15) 
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so that 


and 


„» _ a 0 + a 4 

s lla ~ 1 + b, 


- C ±1U 


’ 21a 


1 + b 4 ' 


b) Connect a perfect load to port 1 and a perfect short to port 2. Then 

s ~ \ o 0 

b - [ 0-1 


so that 


fln Q i 

m _ _ 0 4 


’ lib 


1 - b 


and 


.m _ c 0 c 4 
s 21b ~ 1 - 5* 


Solving simultaneously again yields 

U O’tla + ~ 1 "lla >in, 


a 4 = 


and 


b 4 = 


m m 
5 lla s llb 


2 _ a _0 (Jjja + s Jlbl 

s Ua ~ s Ub 
m 


with similar expressions for c ^ and b ^ in terms of (s^ fl ) and (s^) 


(16a) 

(16b) 

(17) 

(18a) 

(18b) 

(19a) 

(19b) 


3.2.4 Step 4 

a) Connect perfect opens to ports 1 and 2. Then 

c - 10 I 

b ~ [01 J 

so that 

_ (a 0 + a l + V + a 5 _ 0 + a 5 
lla J1 +6^ TTp + b j [i + b ^ 


( 20 ) 


(21a) 


18 


and 

m _ ( *0 + c l + c jl + c 5 

$21 a J1 + b ^ + 5^7 "b 

Note that a, ft and 7 are known. 


7 + «r 

/rr^- 


b) Connect a perfect short in place of the open at port 2. Then 

c - / O' 

~ 0 - 1 

so that 

m _ ( a 0 + a l ~ V ~ a 5 _ a ' ~ a 5 

s ub ~ (7 + - r 5 - w - r 5 

and 

» _ Jo + C 1 - V - c j _ T- - Cf 

s 21b ~ J7 + 1>i - VJ) - 6 ^ ft' ~ 


Eqs. (21a) and (23a) can be solved simultaneously for and 6^ as 


a 5 = 


Ala + « An - ^ + A Ala Alt 
Ala ~ Alb 


and 


a + a- - 0 s n Ua - 0' s 


*« = 


m 

lib 


s lla $llb 


Similar equations for and b ^ result from solving (21b) and (23b) 
simultaneously. 


(21b) 


( 22 ) 


(23a) 


(23b) 


(24a) 


(24b) 


3.2.5 Step 5 

a) Connect ports 1 and 2 with a perfect through line having zero phase shift. Then 


so that 


0 1 
1 0 



a Q + (“2 + a jl ~ a 5 

1 + (^2 A- b 3) ~ b ^ 


q t (^2 
ft + W 2 + b s ) 


(25) 


(26a) 
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and 


so that 


and 


si 
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+ (Cg 
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r + (c 2 
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(26b) 
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~ b 5 

V + (bg 

+ bg) • 

ports 1 

and 2 

with a 
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through line have 180° phase shift. 

Then 



' - [ 
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-1 0 \ 
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~ b 5 " 

V 

+ 

■o 

1 

^Zoo, ) 


c 0 

~ (Cg 

+ Cg) - 

- 

r - 

- (Cg 

+ Cg) 


s 21b 

~ T 

- (bg 

+ bg) - 


P' - 


+ bg) • 

^ZoU ) 


Solving (26a) and (28a) simultaneously for (ag + a g) and (bg + bg) gives 

_ »' °Ua + • s Ub - <r + n s na s "m 


a 2 + *3 - 


and 


m m 
s lla s llb 


a + o' - fi s m Ua - 0' s 


b 2 + b 3 


'lib 


m w 

s 11a s llb 


(29a) 


(29b) 


Similar equations for (eg + Cg) and (bg + b g) result from solving (26b) and (29b) 
simultaneously. Resolving a 2 and a 3 , b 2 and bg, and c 2 and c 3 from their 
respective sums requires some computation, and is discussed next. 


3.3 Solution for a g, a 3 , bg, b 3 , Cg, and Cg 

The calibration measurements of the previous section yield values for the 
following coefficients and sums: a Q , ( a 2 + a 3 ), a ^ a 5 , b ^ (bg + bg), b^, 

bg, Cq, Cj, ( Cg + c 3 ), Cj, and c g. The goal of this section is to express Og, a 3 , 
bg, bg, Cg, and Cg in terms of these known quantities. 
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We begin by returning to the definitions of the o ? , b ■, and c . in terms of the 
network analyzer model [Eqns. (6) and (7)]. Since each of the coefficients of (6) 
may be divided by either s gj or to obtain the form of (7), we will simply let 
Sgj and s s ?8 be given by unity in (6). In effect, we are normalizing to Sgj (or s^) 
for every path leading away from port 1 (or port 8) of the network analyzer model. 

For convenience, the coefficient definitions are repeated here, incorporating the 
normalization discussed above: 


a 0 

= 

4 





(30a) 

a l 

= 

4 4 

~ S 31 

s 44 



(30b) 

a 2 

= 

s s 
s 34 S 51 

~ 3 S 1 

s 54 



(30c) 

a 3 

= 

s s 
s 35 S 41 

~ 41 

<5 



(30d) 

°4 

= 

s s 
3 35 3 51 

~ S 31 




(30e) 

a 5 

= 

s s 
s 3i s 44 

5 55 “ 

5 55 s 34 

s s s 

s 4 i 3 44 s 5i 

s 35 



- 

s s 
S 31 s 54 

s 45 + 

s s 
3 45 s 34 

S 51 + S 41 s 54 

s 35 

(30Q 

b 0 

= 

1 





(31a) 

b l 

= 

3 24 3 41 

" s 44 




(31b) 

b 2 


s s 
s 24 S 51 

~ s 54 




(31c) 

b 3 

— 

S S 

3 25 3 41 

~ s 45 




(31d) 

h 4 

— 

s 25 S 51 

-4 




(31e) 

b 5 

= 

s s 
3 44 s 55 

-4 

S S 

s 24 S 41 

S S S 

s 44 S 51 3 25 




- 

5 S 

3 54 45 

+ s 45 

s 24 S 51 

. s s 

+ 3 41 3 54 25 


(310 

c 0 

= 

3 S 

S 61 





(32a) 

C 1 

= 

s s 

s 64 41 

- S 61 

4 



(32b) 

c 2 

= 

s s 
s 64 51 

- S 61 

4 



(32c) 

c 3 

— 

S S 

s 65 S 41 

_ 

S 61 

4 



(32d) 

c 4 

= 

s s 
3 65 3 51 

_ 

S 61 

4s 



(32e) 
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and 


c 5 ~ 


S 61 S U S 55 ~ s 55 s 64 S 41 " s 44 S 51 s 65 


~ S 61 s 54 s 45 + s 45 s 64 S 51 + S 41 s 54 s 65- 


(32t) 


At first it may not appear to be possible to solve for the subscripted constants s • 

U 

separately. However, the complete solution may be considered a possibility since we 
have already obtained 14 measured quantities on which the 14 coefficients depend. If 
we look even further, we see that we may also set to unity if desired since it 
always appears in product with other terms and may be absorbed. 

To consider the solution for the additional coefficients, we begin by combining 
Eqs. (30) - (32) to obtain 

(33a) 
(33b) 


a 2 a 3 ~ 


b 2 b 3 ~ 


a l a 4 a 0 a 5 
b l b 4 ~ b 5 


and 


with 


c 2 C S ~ 


C 1 c 4 c 0 c 5 


°2 + a 3 ~ 

b 2 + b 3 = 


'a 

u 


(33c) 

(34a) 

(34b) 


and 


c 2 + c 3 = V (34c) 

Since and (and similarly, a^, 6^, a g, and by) are the only unknowns in 

(33c) and (34c), we may solve for them as quadratic forms as 

7, 

(35a) 


c t = 


~T * ~T + ~^4 “ * c o e t 


and 


c 3 = ~T ~ ( * ~r R + 4c t c 4 - ~4 c^ 5 ). 
Similar forms occur for a^, a ^ b 2 > and b ^ 


(35b) 
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The resolution of the (±) is based on the application of physical requirements 
for the system, specifically that is dominated by the actual s,^ of the calibration 
standard. To this end we may equate with the larger of the roots in (35) and 


with the smaller. The remaining coefficients require a bit more work, 
by rewriting combinations of (30) through (32) as 

We begin 

a l ~ b l a 0 = S 4I ( s 34 ~ 8 24 a 0 ) 

(36a) 

a 4 ~ b 4 a 0 = S 51 ( s 35 ~ s 25 a 0> 

(36b) 

C 1 ~ b l c 0 = s 4l ( s 64 ~ s 24 c 0 ) 

(36c) 

and 


c 4 ~ b 4 c 0 = S 51 ( s 65 ~ s 25 c 0 } ' 

(36d) 

We may also combine (34a-c) to obtain 



7 a 7 b a 0 ” S 51 s 24 a 0 ') + S 41 ( s 35 s 25 a 0 ) 

and 


7 c ~ 7 b a 0 ~ S 51 ( s 64 s 24 c 0 ) + S 41 ( s 65 s 25 c 0 ) 

into which we may substitute (36) to obtain 



7 a - 

~ 7 b a 0 

S 51 ( 
S 41 1 

‘ b l a 0> + 

*JL 

s 5l 

( a 4 b^a Q ) 


(37a) 


7 c ' 

- 7 b c 0 

= 

s 4l 1 

+ 

S J1 

S 51 

( c 4 ~ b 4 c 0^‘ 


(37b) 

Eqs. 

(37) provide two forms for s ki/ s ij 

given by 




S 51 

( y a 

- n a o><u - 

h fo> - 

( y e ~ 

■ h c o>( a 4 - 

b 4 a o) 

f38al 


s 4i 

7^7 


w ■ 

(°1 ~ 

■ b i c o)( a j - 

b 4 a 0 j> 

l UOd 1 

and 










S Jl 

(“1 

1 

<3 

1 

b 4 c t>) " 

( C 1 ~ 

■ - 

b 4 a o) 



s 4i 

7* 7 


w - 

( C 1 “ 

' b l c 0>(la ~ 

7 b a o)’ 

(dOU J 
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If we now substitute (32a and b) into 32(c) we obtain 

c 2 = C 1 + c 0 (jjj s 44 ~ s 54 )* ( 39 ) 

Only the last term of this expression is unknown at this stage and we may write it 
as 


<*;; 3 u v - 


and similarly for the term 


(i il s 

{ *51 55 


"V = 


s _51 C 1 
s 4l c 0 


til -Hi, 

s 5l c 0 


(40a) 


(40b) 


Forms for the unknown coefficients a^, a bg, and by similar to (39) may be 
written and (40a and b) used to obtain 


— ( 2J. \ n i a _0 i . S 51 \ 

^s,) a l + c Q ( c 2 s 41 c l ) 


*51 

Ui 

= (i^) a 4 + 

S 51 4 

= (t 1 ) bj + 
s 4i 1 


d mm 

c 0 ' 3 S 51 C ^ 


o _i a 


i 


51 


c Q 




s 4l ^ 


and 


= (j^) h + 

s Sl 4 


( c s ~ s*r C J- 


(4U) 

(41b) 

(41c) 

(4 Id) 


0 " 51 

This completes the solution for the coefficients of (30) to (32). An analogous 
procedure may be used for the coefficients with an input at port 8. 


3.4 Calibration Using Imperfect Standards 

The ideal standards envisioned in Section 3.2 do not exist in practice. To the 
extent that the actual standards differ from the ideal, the task of calibrating the 
network analyzer becomes more difficult. 
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The simplest non-ideal case occurs when some or all of the non-zero entries in 
the S matrix differ from their ideal values of +1 or -1. An example is a standard 
whose electrical length is incorrect, resulting in a phase angle other than 0° or 
±180°. Mathematically, this case causes little difficulty: the equations relating the 

measurements [( 5 ;/ fl ) an( I (tfu,)) or [( 5 2/ a ) an( * ( s ^ )] to the coefficients a- 
and b- or c and b ■ generalize to the form 

s a = § <«a) 


and 


s b 


a' + 7 'x 
V + d'p 


(42b) 


where a, 0, 7, 6 , a', /?', 7', and S' are known and x and y are unknown. These 
equations have the simultaneous solution 


x 


a' $„ 
a 


4~ s _b - (0' --T&) \ s _b 

- r * fl + -f 7 ^ 


(43a) 


-4r [~ 0 7' + <*' 7 + 0 Y s - 0' 7 sj 
y = ^ - (43b) 

- r s a + ~r 7 s b 

It is easy to check that the equations of Steps 1 through 5 in Section 3.2 are special 


cases of (42a)-(43b). 

For example, in 

Step 2a 

we can 

identify 

s = s m 

s a s lla’ 

a = a Q , 

7 = 

1, 

x = 


0 = 

s = 

1, 

y = 

and in Step 2b 





, - o m 

3 b ~ s llb’ 

a' = a Q , 

7 ' = 

-1, 

x - 


0' = 1, 

S' = 

-1, 

y = 


25 


The next level of difficulty occurs if the formerly zero entries of S are allowed 
to be nonzero. One example is the use of an imperfect load having a nonzero 
reflection coefficient; another example is a through line with discontinuities that 
make 5^ or 5^ nonzero. In this case, the procedures of Section 3.2 are no longer 
applicable. Instead, the nine test conditions 1 through 5b yield nine simultaneous 
equations in a - and b it and another nine equations in c - and b •. Solution for the 
unknown a- and 6-, or for the unknown c- and 6 -, requires solution of a ninth 

XX XX 

order system of equations. Still, apart from the extra computation there is little 
difficulty in dealing with this case. The underlying (and perhaps dubious) 
assumption is, of course, that the S-parameters of the standards are accurately 
known. 

The highest level of difficulty comes about when the characteristics of some or 
all of the standards are not completely known. An open, for example, might 
reasonably be expected to have a reflection coefficient of magnitude one (implying no 
loss), but its phase angle might well be nonzero due to fringing of the electric field 
at the end of the transmission line. Our proposed techniques for dealing with this 
last situation are described in the next chapter. 
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CHAPTER 4 
FULL CALIBRATION 


The previous chapter discussed ways of calibrating the network analyzer when 
known standards are available, whether or not those standards are "perfect." This 
chapter takes up a more difficult problem — how to calibrate the system when the 
characteristics of the standards are not completely known. We wish to emphasize at 
the outset that this statement does not imply that the network analyzer can be 
calibrated accurately using inferior or damaged standard components. Certain 
assumptions based on their physical properties must be made in order for our 
proposed method to work. Nevertheless, this method represents a significant 
theoretical advance over that previously recommended by Hewlett-Packard. 

4.1 Standards 

We begin by discussing the standards that will be required for this calibration 
procedure. These standards are all sold by Hewlett-Packard for calibration and 
verification of their network analyzers. Some of them are not supplied as part of the 
"basic" calibration kits, however. Thus, implementation of our recommendations 
might require purchase of additional components that are, admittedly, not 
inexpensive. 

Fortunately, the full procedure described here is necessary only when the 
standards are not completely known. Once the standards are completely 
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characterized, the abbreviated calibration procedure of Chapter 3 is adequate provided 
the standards do not change in their characteristics. 

The process can be carried one step further (in face we recommend doing so) 
by using the full calibration procedure and a set of prime standards to characterize a 
second set of working standards, to be used for day-to-day calibration with the 
abbreviated procedure. Thus, one set of prime standards can be maintained for an 
entire organization employing many network analyzers. An additional advantage is 
that the prime standards will receive infrequent use, hence will tend to remain in 
much better physical condition. 

To illustrate our calibration scheme, we assume for the remainder of this 
chapter that the goal is to characterize a set of prime reflection standards and in the 
process calibrate the network analyzer for reflection measurements at port 1. Using 
reflection measurements, the working standards can then be characterized. 

The required prime standards are as follows: 

1. Perfect short: As before, we specify that the component have reflection 

coefficient r = -1 + jO. The phase of 180° is an arbitrary choice to provide 
a phase reference for the system; in effect this choice locates the reference 
plane precisely at the shorted end of the transmission line. Since there is no 
fringing of the TEM wave on a transmission line in the vicinity of a short, 
there should be no frequency-dependent behavior. The unity magnitude of the 
reflection coefficient is based on an assumption of losslessness in the component 
(including its connector). A worn, damaged, or dirty connector will render this 
assumption invalid. Further, there is no direct way to test the short, since it 
is being used to calibrate the network analyzer in the first place. For this 
reason, Hewlett-Packard recommends that all of its precision microwave 
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standards be verified using sophisticated physical, rather than electrical 
measurements. 

Since some loss, however small, is inevitable in a real-world connector, 
there may be a legitimate concern about the lossless assumption. The best 
justification for this assumption is that no better one can be made unless the 
short is modeled with the material properties used in a skin-depth formulation. 
This modeling would be based on documented estimates of the material 
conductivity obtained from cavity resonance measurements. Further, if the 
physical quality of the standards is very high, it seems reasonable to assume 
that the loss associated with each mated connector pair is about the same. 
Under this assumption, the loss of one connector pair can be absorbed into the 
network analyzer model, so that it is the component beyond the connector that 
is being characterized. 

2. Lossless open: We assume that the open has a reflection coefficient close to 

unity which is consistent with the existing closed region. Because of field 

fringing near the discontinuity, the precise electrical location of the open end is 
not known in advance and it exhibits a frequency-dependent behavior. As with 
the short, the unity magnitude of the reflection coefficient is a consequence of 
its assumed losslessness. 

3. Two lossless transmission lines of different electrical lengths: Once again, the 

lossless property is assumed. There is no need to know the exact lengths of 
the lines; it is only necessary that they be significantly different. 

4. Sliding load: This component consists of a transmission line terminated by a 

block of resistive material. The termination may be slid up and down the 
transmission line to change its distance from the input connector. Assuming 
that the termination is of good quality and that it maintains good contact 
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with the transmission line, the reflection coefficient seen at the input of the line 
will have a constant magnitude. Its phase, however, will vary with the position 
of the terminator (load). Thus the input reflection coefficient of the line will 
be located on a circle in the complex plane. 


4.2 Reflection Calibration 

Using the precision standards just described, the network analyzer is to be 
calibrated for reflection (sjj) measurements. We assume that port 2 is terminated 

with a good matched load so that s ^ = s ^ — 0 (Strictly speaking, this 

termination is not a requirement, as non-zero s^* Sgp and S 22 may be absorbed in 
the constants as long as they do not change.) The describing equation of interest is 
(6a), which reduces to 
m 


a o + a i s n 
1 1 7 + 


(44) 


4.2.1. Shorted line measurements 

With the short connected to port 1, 


and 


’ll 


= -1 


s'" 

5 11 


°0~ a l 

t~^t: 


= r 


(45). 


If the short is now connected to port 1 through either of the two lossless 
transmission lines, the input of the transmission line must have s.j of the form 

s n = 1 JV, (46) 

that is, a reflection coefficient of unity magnitude (because of the lossless nature of 
the line) but nonzero phase shifts (because of the nonzero lengths of the lines). 
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Obviously, the short alone is a special case of (46) having 


The measured s 


tn 

11 


short 


S 11 


short 1 


5® 

5 U 


short2 


« = a 0 + V'*’ 1 ' 
111 1 + 


<p = 0. 

satisfy the equation 


47) 


for each i. It is possible to show that the bilinear transformation relating to 
Sji maps circles into circles; specifically, if 

s n = |T| e jd , 0 < 0 < 2j (48) 

then 

s m n = C + R e ja (49) 

where 


1 - 


C = 


l l 


^ 1 + 1 


^7 


bj \ 2 |r | 2 


(50a) 


R = 


a 0 b l . 


(50b) 


1 - i^i 2 |T| 2 

Thus, connecting the short through various lengths of lossless line produces 
measurement points that lie on the circle given by (49)-(50) with |T| = 1. A 
minimum of three points is needed to determine the circle uniquely. These points, 
as discussed above, can be found by measuring the short alone, and the short with 
two different lengths of lossless line. (As an alternative and improvement, the short, 
the open, and both the short and the open on a lossless line can be used for 
determining the above circle. This approach circumvents some of the problems which 
arise due to loss in the line extensions.) 

The circle in the s^ plane that results from these measurements will be 
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denoted 


with 


and 



C+R 
s s 



n 1 

aj a 0~~5 T 
s ~ ~ r i + 1 - 16^2 

a _0 ~ a _l b _l 
1 - \bj\ 2 


A l 


b _J_ a jl 

1 ~ \bj\ 2 ' 


4.2.2. Sliding load measurements 

At the input of the sliding load 

s" n = e e™t 


( 51 ) 


(52a) 


(52b) 


(53) 


where 0 < e < 1. Measuring the sliding load for various positions of the 
termination gives values of satisfying 


Qq + aj e e^l 
1 + bj e 


(54) 


for all ^'s. Applying (48) - (50) gives the center and radius of the measurement 
circle in the plane as 

s "iu = c t + R t e3a t (55) 

with 


a _0 ~ a J b j_ ^ 

1 - |6 ; | 2 e 2 


(56a) 
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and 


\a Q bj - a t \ e 

1 - l^l 2 < 2 * 


(56b) 


Again, at least three measurements (three positions of the termination) are required 
to determine and R p. 


4.2.3. Solution for unknown coefficients 

Eqs. (45), (52), and (56) involve the unknown coefficients dg, a p and b ^ the 
unknown sliding load reflection magnitude e, and the known measurement quantities 
r , C , R , C f , and Rp. We note that there are eight real equations — two each 
for the complex-valued equations involving T,, C ,, and C», and two real-valued 
equations involving R $ and R^. Since there are only seven real unknowns (one for t 
and two each for dg, a^, and bj), one equation is redundant. The equations used 
for the solution enforce a) the short to be on the unity circle (an assumed lossless 
condition), b) the computed reactance circle to be a circle of unity radius centered at 
zero, and c) the load circle to be centered at zero. These equations may be solved 
exactly to obtain 


dg - C, C s / 2 + 


e v s 


* i -*} 


and 


- \((t] - «}) 2 + (\C(-C S | 2 + Up) \c r c s I 2 ; 

/ tnCf-c/], 

= <- c o - “ o > (r, - V / K + < r s - 9 V> 

= - (*J - C o <«. - (b - «0> 

/ <3 + < r s - c 0> c V ««)) 


(57a) 

(57b) 


(57c) 
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e = (o g bj - aj R ( / {{bj C ( - aj 2 - b] Rj). (58) 

4.2.4. Determining radius and center for the reflection circles 

The short and sliding load measurements of the previous sections determine two 
circles in the complex s ^ plane, 

5 7is ~ C s + R s e} ° s 
and 

5 T a = c i + K i eJat 

respectively. Here we investigate the mathematics of computing the radii and 
centers. 

Suppose first that we have three measurements Sp s and s^, and wish to 
compute R and C for the circle on which they lie. We begin by drawing two chords 
joining these points, say chord L ^ joining Sj and and chord joining and 

s^. The midpoints of these chords, say p ^ and p^ respectively, are 


P 1S - 


S 1 + s t 


Pis - 


S 1 + S s 


The perpendicular bisectors of L ^ and say B ^ and 21^ respectively, pass 
through and and intersect at the center. We thus need to find the 
equations of B ^ and to determine where they intersect. These bisectors are given 
by 

B 12 = p 12 + °12 tt s l ~ s 2 ) ( fi 0a) 


B 23 ~ p 23 + a 23 J ( s 2 s 3 ) 


(60b) 


where a 19 and a 9<) are real numbers. If we equate B 19 and B 9q , the a’s may be 
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determined to give the center C as 


C = 


[s 1 (\s 3 \ 2 -\s 2 \ 2 ) + s 2 (\s } \ 2 -\s 3 \ 2 ) + Sgds^-lsjl 2 )] 
l 3 I s 2 ~ s l s 2 + S 2 S 3 ~ S 2 S 3 + s 3 s 1 ~ S 3 S 1 3 


(61) 


The radius is then given by | - C |. The best numerical conditioning results when 

the three points Sp anc * are nearly equally-spaced around the circumference 
of the circle. This conclusion is true for two reasons: 1) with (s^-s^) or 

small, the denominator of (61) becomes small, leading to possible roundoff error. 2) 
with (Sj-Sg) or (s^-s^) small, R and C become very sensitive to small errors in the 
measurements (due, for example, to noise in the system). 

If the measurements s • are noisy, or are suspect for any reason such as lossy 
connectors, it is preferable to use more than three points to define the circles. In 
this case, a best-fit procedure, rather than Eq. (61), is needed. One possible best-fit 
approach is described in the next section. 


4.2.5. Least-squares circle fitting 

Our problem is to find the circle that "best" fits a collection of measurements 
{s • | 1 < i < k, k > 5), in some sense of the word "best." A common measure of 
the goodness of fit, and the one that will be employed here, is least squares (the 
curve, in this case the circle, is chosen to minimize the mean square error in radius 
over the data points). 

Once again, the measured points in the complex plane will be denoted 
Sp s n • The desired representation is a circle 

s u = C + R e- Jlp . (62) 

The error between measured point s ■ and the best-fit circle is 

U = \(\ Si - C\) - R\. (63) 
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Possible phase errors in s- (that is, errors that displace s ■ around the circle but not 

away from it) do not influence the best-fit choice of C and R. Thus the 

least-squares problem reduces to that of choosing C and R for minimum error. 

Using the error representation of (63), the mean square error is given by 

F(C,R ) = mean square error 

N 

= l (|«,. - Cl -If. (64) 

i=l 

To determine the minimum of F, we differentiate with respect to R, Re(C), and 

Im{C) and set each derivative to zero. Thus we have 

N 

JR ~ ~ 2 I d 5 *' - ^1 “ *) = 0 (65a) 

i=l 

dr v , , , R *( c - s i) 

MefT) - 2 l ~ c \ ~ | e = 0 ( 65l >) 

i=l 1 

and 

dF v Im ( c ~ s j) 

JTm{UJ = ~ 2 l(\ s i ~ °\ ~ R ) - |7T- 5 . | = °- ( 65c ) 

Eq. (65a) is easily solved for ^ to obtain 

N 

R — -j- \ l 5 z - ~ £| > (66a) 

i=l 

while Eqs. (65b) and (65c) may be combined to obtain the function 

N C-s. 

/ = C -r -p- £ (s • + R |£ _ /|) = 0. (66b) 

i=l 1 

In general, (66) must be solved iteratively. We have implemented a secant 
optimization routine which performs well for most cases of interest. On occasions, 
the center will tend to infinity, which also offers a zero to f. This tendency is 
associated with clustering of the data points with a related inability to. fit a circle 
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properly. This problem specifically occurs at 10.7 GHz for the sliding load used in 
the HP system with the spacings designated on the sliding load body. Even HP 
recommends additional sliding load data points to alleviate this problem. 

To help avoid the potential infinities, F was modified to include an additional 
R 2 multiplier as 

N 

r = l (|s, - C\ - R) 2 R 2 . (67) 

i=l 

To simplify the process, R was taken to satisfy (G6a) rather than the new equation 
obtained by differentiating F with respect to R. This choice gives an optimizing 
function of the form 

N C _ S ' 

f = [C “ 7 - J («,• + * [77 _ 5 *|)7 & 

i=l 1 

N N C-s. 

+ 2 R ~T a £ (I 5 * ~ ~ *) 2 I (|<7 _ a *|)= °- ( 68 ) 

i=l i=l 1 

Again, (68) is solved iteratively for C using (66a) to define R. This equation has 
resulted in excellent results for the computed circles. Hewlett-Packard has suggested 
a modified procedure for finding the circles which minimizes the mean-squared error 
of the squared radius. A definite advantage of this procedure is the existence of a 
closed-form solution requiring no iteration. This approach has been modified to 
include the reactance circle form to be presented in Sec. 4.2.6. The resulting 
equations have been implemented for the calibration process and has been found to 
give excellent results. An alternate form has been developed which minimizes the 
least square error of the center from the bisectors of all point sets. This latter 
technique may be solved in closed form and gives excellent results as shown in Fig. 4 
for the phase of the HP "open". The figure also includes the 
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Figure 4. a) Calibration of the HP "open" phase based on the modified HP 
circle fitting procedure and the bisector form, b) The effects of the 
sliding load and a fixed load on varied line lengths are also shown. 


38 


effects of using the sliding load versus fixed terminations mounted on the ends of 
airlines, to be discussed in more detail in the next chapter. The modified HP form 
is currently being used for the curve-fit procedure. 

4.2.6 Center and radius of reactance circles 

The reactance circles provide for a modification in the circle fitting procedure. 
First, we assume that the short is ideal and the circle must go through this point. 
As a result it is simple to estimate the radius as 

R = | s t -C\. (69) 

In addition, we assume that the "open" is an ideal reactance. In other words, it has 
a reflection coefficient of unity magnitude. Thus the center of the reactance circle 
must lie on the curve defined by 

C = + j x (s jr — 5^) (70) 

with x the only unknown. 

If we now differentiate (67), after the substitution of (69) and (70), with 
respect to x and set the result to zero we obtain 

/ = -JT = * * 2 1 1 (/ “ Tv=tt) Ser{C ' s ”) U {s ‘ " Sg)),] 

+ [{ 1 - B j T 2 + (1 - t^B J ^ )] |s; - i^| 2 zj = 0. (71) 

With this form, the algorithm uses a short circuit reference, and a collection of 
additional reactance references (not necessarily of known reactance). The open circuit 
is used as a preferred first unknown reactance on which the line defining the circle 
center is based. This is primarily used due to the probability of additional loss in 
the other references. The remaining reactances are obtained by using several airlines, 


39 


preferably with nonharmonic related lengths, terminated in both a short circuit and 
an open circuit. The used of both the short and open is to counteract the errors in 
the circle determination arising due to loss in the airlines. The short and open 

combinations tend to average out the loss effect by appearing on opposite sides of 
the reflection coefficient plane. As with the load curve-fit, the reactance circle 
routine has been modified to use both the bisector and the HP approaches rather 
than the minimization of (67). These approaches have given improved results and 
the modified HP form is currently in use. 

Detailed discussions of the calibration process and the measurements of 
component standards are provided in the next chapter. In addition, models are 
developed for the actual open circuit references and terminations. Consideration is 
also given to the possible imperfections in the short circuit reference. 
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CHAPTER 5 

RESULTS OF REFLECTION MEASUREMENTS 


This chapter presents the results of the calibration procedures developed in the 
previous chapter for reflection measurement. This procedure requires no source or 
test port switching and thus does not have any errors associated with such switching. 
The errors which occur are thus fundamental errors in calibration or system errors 
associated with the network analyzer sweep circuits, frequency sweep, or 
discretization. The data was taken on an HP 8510 network analyzer system with an 
Hp 85080A, 7mm, calibration kit. Initially, we will consider the basic measurement 
and associated calibration. We will then consider the modeling of both the "open" 
standard and the fixed termination. 

5.1 Basic Measurement and Calibration 

In the IIP system, two termination standards are provided for use. An 
excellent fixed termination is provided for use below 2 GHz and a sliding load is 
provided for use above 2 GHz. To avoid problems with imperfections in the fixed 
termination, it was used at the end of a series of 50 ohm, beaded, precision airlines 
to provide a load circle similar to the sliding load. This approach to the 
termination problem gives good results over the entire frequency range as shown in 
Figure 5. However, the fixed load does have an oscillatory error. In the realm 

above 2.0 GHz, the sliding load provides excellent results. In fact, the sliding load 
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Figure 5. Comparison of the calibration for the HP open phase using a 

fixed termination terminating airlines and using a sliding load. 
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works well down to 0.5 GHz. The reactance calibration was obtained with both 
short and open references and both the short and open terminating the same airlines. 
The open was treated only as an unknown reactance while the short was considered 
perfect. The data calibration is an improvement over the simple calibration using 
the termination references, the HP open, and a short. The full calibration required 
the determination of circles which best fit the calibration data discussed in 
Chapter 4. 

5.2 Modeling the Open Calibration Standard 

To address the full calibration process, we are interested in determining a 
model for the HP open. We have three different models for the open as a frequency 
dependent capacitor readily available. A modified form of the 2— term data for the 
lip open #2 is provided by Innovative Software. This model is a two term power 
series given by 

C = 0.0 905 + 0. 0000785 p pF (72) 

for / in GHz, with a phase error for the resultant phase calibration as shown in 
Figure 6. As in all the plots • to follow, the calibration has been done using the 
short reference, airlines terminated in both a short and a open for the reactance 
circles, and either the sliding load or a fixed termination on a series of airlines for 
the 50 ohm reference. The estimated phase in Figure 6 was determined by using the 
C of (72) to model the open. 

Similar results are shown in Figures 7 and 8 for three and four term power 
series approximations to the capacitance given by Hewlett-Packard. These series are 
respectively 

C = 0.09285 * 0.0000072 p * 0.0000043 p pF (73) 
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Figure 6. Comparison of the calibration for the HP open phase to the 

phase based on the two-term capacitance model of IS. 
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Comparison of the calibration for the IIP open phase to the 
phase based on the three-term capacitance model of IIP. 
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Figure 8. Comparison of the calibration for the HP open phase to the 

phase based on the 4-term capacitance model of HP. 
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and 


C = 0.0872 * 0.001695 f - 0.00015 p + 0.0000089 p pF . 

( 74 ) 

The initial goal of this work has been to develop a model for the capacitor. 
This improved model may then be used for future calibration. The process used to 
calibrate the open may also be used in situations when the open may not be used 
directly. To model the open, the calibrated phase (or effective capacitance) was fit 
to a corresponding model of the open. This curve fit was forced to minimize the 
phase errors between the measured data and the model. Four basic models were 
chosen for consideration: a two term series like (72); a three term series like (73); 

a three term series (denoted formula) with the series forced to the capacitance of the 
HP three term scries at 40 GITz; and an inductor in scries with a capacitor to 

represent the delay inherent in the open circuit connection. The results are 
indistinguishable in plot form as shown in Figure 9 for the inductive expansion and 
the first 3— term series. 

To obtain the results of Figure 9, the phase of the open reference was 
matched to that of the model used for the open. The results for the inductive 
termination model are shown in Figure 10. This model shows an improvement over 
the fits of the HP 3— term model and the Innovative Solutions model shown in 

Figures 6—7. A similar improvement is also found if the calibrated capacitance of 
the open is compared to the estimates based on the other two circuit models. 

Figures 11 and 12 depict the calibrated capacitance using the reactance circle 
approach as compared to the HP 3— term model and the inductive model. In 

particular, the capacitance is over estimated for the HP model. The low frequency 
error arises due to the inability to fit a circle to the load points adequately. 
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Figure 10. Comparison of the calibration of the HP open phase to the phase 
based on the Virginia Tech inductive model. 
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Figure 11. Comparison of the calibrated capacitance of the HP open and the 
HP 3— term series model. 
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Figure 12. Comparison of the calibrated capacitance of the HP open and the 
Virginia Tech inductive model. 
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The importance of the calibration standard modeling is emphasized in a plot of 
the reflection coefficient of a 30 cm shorted airline. This is the classic problem 
which has made researchers suspicious of the HP calibration system. If the 
standard models arc in error, the response for the 30 cm airline shows an oscillatory 
behavior which may even give the false concept of an airline with power gain. Such 
data is presented in Figures 13 through 17 using data collected on an HP 8510 
computer-based network analyzer system.* The HP 3— term and 4— term series are out 
of phase with those of the Innovative Software data at the high frequency end and 
represents an apparent underestimating of the capacitance of the open standard for 
the Innovative Software data. The Virginia Tech calibrations using a 3-term series 
and an inductive model both fit the data with large excursions. Though there was 
and excellent fit to the open data, there appears to be an error in the Virginia Tech 
capacitance estimate of the open. However, this data was taken for a 30 cm beaded 
airline which may have some resonances due to the beads. To check our results, we 
should view the 30 cm line terminated in a short circuit using the sliding load data 
of Fig. 18. This data used the collection of airlines for calibration rather than the 
open model. This data is compared with the 3— term VA Tech model and the 
4— term HP model in Figs. 19 and 20 respectively. 

5.3 Modeling the Fixed Termination Standard 

In many instances, we would expect the calibration effects of the termination 
reference also to become important. To develop a model for the termination, we 
may consider the amplitude and phase plots of Figure 21. Though there is a fair 
amount of noise on the data due to the resolution which is approaching the 
discretization of the system used, it is apparent that the termination may be modeled 
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Figure 13. Calibration of a 30 cm beaded airline using the HP 4— term 
capacitance model. 
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Figure 15. Calibration of a 30 cm beaded airline using the Innovative 
Software capacitance model. 
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Figure 16. Calibration of a 30 cm beaded airline using the Virginia 
Tech 3-term capacitance model. 
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Figure 17. Calibration of a 30 cm beaded airline using the Virginia Tech 
inductive model. 
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Figure 19. Calibration of a 30 cm beaded airline on a HP 8510 using the 
VT 3— term capacitance model and the airline technique. 
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Figure 20. Calibration of a 30 cm beadless airline on a HP 8510 using the 
HP 4-term capacitance model and the airline technique. 
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Figure 21. Calibration of the fixed termination using the inductive model for 
the open and the fixed termination with several airlines to define 
a load circle. 
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reasonably well by a fixed resistor at the end of a length of line. The actual length 
is difficult to estimate, though it appears to be between 3 to 4 cm. A description of 
the actual termination is shown in Figure 22, with the appearance of a thin film 
resistor 0.8785 inches from the connector face. Using a HP 4 192 A low frequency 
impedance analyzer at 1000 Hz, the resistance is found to be 49.30 ohms at -0.12°. 
This measurement gives a reflection coefficient of -0.007 which differs from the 

c 

—0.002 reflection coefficient below 2 GHz of Fig. 22a, though within reasonable limits 
in absolute resistance. This latter corresponds to 49.9 ohms at low frequencies. If 
we estimate the length to the resistance to be the 0.8785 in distance to the 

connector face, the one half to account for its distributed nature, we obtain an 
effective length of 2.231 cm. The data for 49.8 ohms located 2.231 cm from, the 

connector face has been plotted in Figure 21. The result is about twice that needed 
to fit the data. The error results from the capacitive washer which effectively 

shortens the line. 

The effect of this termination model is not always critical as is seen in the 
reflection amplitude of the open standard. These results shown in Fig. 23 were 

obtained assuming the termination was a perfect 50 ohm termination. The 
calibration of the open circuit reference using the 50 ohm assumption has excellent 
results below about 5 GHz. 
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Figure 22. Construction details of the HP fixed termination standard. 
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Figure 23. Calibration of the 30 cm beaded airline using a) the fixed 
termination on a series of airlines and b) the fixed termination 
exclusively as a 50 ohm standard below 2 GHz. 
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5.4 The Effect of a Nonideal Short Circuit Reference 

Before we complete this chapter, we should consider the last calibration 
assumption used in commercial systems in addition to the calibration steps above. If 
you have not already considered the possibility, the short circuit standard does not 
necessarily have a reflection magnitude of unity. We now address this last concern, 
but we will not complete the implementation necessary to incorporate the nonideal 
short effect into the calibration in this report. 

We assume that the electromagnetic fields in the vicinity of the short circuit 
reference form a TEM mode in the transmission line. At the short we will 
determine the relection in a planewave sense from a conducting wall with a finite 
conductivity. In the planewave problem, the characteristic impedance of copper is 



For copper p is x 10 7 Il/m, t is approximately x 10‘ 9 F/m, and a is 

5.8 x 10 7 S/m. Determining the reflection coefficient of 7/ c with respect to ?/„ for 
free space we obtain 



which at 10 GHz is ( 0.99986 / 179.99° ). At 40 GHz this becomes 
( 0.99945 /179.97 0 ). Thus we find that the effect of the loss in a copper short 
circuit standard should only cause a variation on the order of 0.05 %. If the system 
is to be taken to the limit for a precision of 0.01%, then this correction would 
become critical. In the applications for which systems are currently used, this 

correction is not necessary and is not considered further here. 

We have seen in this chapter, the importance of proper modeling of the 
standards used in the calibration of a network analyzer system. If such corrections 
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are not used, the resultant data can not attain the accuracy limits of the system in 
use. This is even true for an HP8410 network analyzer system under computer 
control. With the increased system stability and smaller discretization error, the 
limits of the HP8510 network analyzer system require even tighter constraints. 



CHAPTER 6 

CONCLUSION AND RECOMMENDATIONS 


This report has presented the fundamental principles and necessary calibration 
procedures for the accurate measurement of S (Scattering) parameters. These 
important parameters require a series of steps to properly calibrate a system used to 
measure the parameters. This is particularly true for measurements made over a 

broad frequency range for which the usefulness of this parameter set in the design 
process is ideal. Measurements may be made from below 1 MHz to the millimeter 
wave spectrum. We have shown the calibration process and standards to be at the 
heart of such a measurement system. 

We have investigated the calibration procedures used by Hewlett-Packard and 
have found the same anomalies in the resultant data as have other researchers. The 
purpose of this report has been to investigate the full calibration process for the 
existing HP measurement system. We begin by reviewing the basic calibration 

process and measurement system. The approach taken by Hewlett-Packard was 

discussed as the methods were developed. A full model for the system was also 

developed. The model was used to develop a full system calibration, including the 
calibration of the HP "open" circuit reference. This approach has the additional 
advantage of not requiring the "open" for the calibration, a fact critical to certain 
measurement environments. In such instances, the calibration must often be done 
using some alternate procedure such as that presented in this report. More 
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importantly, the basic formalism of the modeling enables a user to develop a 
calibration procedure tailored to a specific need where the HP procedure is 
inadequate. Our premise has been to improve the IIP calibration data and then use 
the fundamental HP procedure when possible. This latter procedure provides a high 
degree of accuracy if the correct calibration data is available, usually in excess of the 
accuracy specifications of the HP 8510 network analyzer system. 

Sample results have been presented which validate the method and 

demonstrate the potential errors using HP's process and data. These results have 
included the effects of both fixed terminations and sliding loads on the choice of 
calibration data used to represent the HP "open" standard. Limitations of the 
approach were discussed in addition to assumptions built into the calibration and 
accuracy process. The final results provide the fundamentals of an accurate 
calibration process for network analyzer systems. Two basic open reference models 
were found to give good results. The first model was a 3— term approximation to the 
open capacitance given by ( 0.09183 + 0.000006 J* + 0.0000046 fi) pF. There should 
be no linear term as used in the newest of the HP models (in fact no odd terms) to 
meet with proper synthesis theory which would indicate such terms would represent 
loss, not capacitance. The 3-order term accounts for neglecting higher order terms. 
The second model consists of a 0.09135 pF capacitor in series with a 0.000205 /d 1 
inductor representing the transmission line effect. These models work well and are 
recommended. 

Additional work is needed to implement of the full calibration process on a 
Hewlett-Packard or other computer system. Also, no actual work has been done in 
implementation of the full 2-port calibration process. An obvious extension of this 
work is to the new measurement environment provided with the HP 8510D network 
analyzer system. This latter system has incorporated the new process mentioned in 
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the introduction which avoids difficulties with the source switching. However, the 
current implementations are known to provide certain fundamental limitations, due to 
imposed assumptions for simplifying the calibration process. These limitations and 
the ultimate effects on the calibration of the data require further investigation. 
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